Critical Behaviour of the Drossel-Schwabl Forest Fire Model 
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We present high statistics Monte Carlo results for the Drossel-Schwabl forest fire model in 2 
dimensions. They extend to much larger lattices (up to 65536 x 65536) than previous simulations 
and reach much closer to the critical point (up to 9 = p/f = 256000). They are incompatible with 
all previous conjectures for the (extrapolated) critical behaviour, although they in general agree 
well with previous simulations wherever they can be directly compared. Instead, they suggest that 
scaling laws observed in previous simulations are spurious, and that the density p of trees in the 
critical state was grossly underestimated. While previous simulations gave p sa 0.408, we conjecture 
that p actually is equal to the critical threshold p c — 0.592 . . . for site percolation in d — 2. This is 
however still far from the densities reachable with present day computers, and we estimate that we 
would need many orders of magnitude higher CPU times and storage capacities to reach the true 
critical behaviour - which might or might not be that of ordinary percolation. 



I. INTRODUCTION 

Empirical analyses suggest that power laws with 
anomalous exponents are ubiquitous in nature, ranging 
from 1/f noise and earth quake distributions to fractal 
coast lines, species extinction rates, weather records, and 
the statistics of DNA ||, |, |, (§• On the other hand, it 
is well known that such scaling laws - most clearly seen 
by linear relationships in log- log plots - can be spurious. 
Log-log plots have a notorious tendency to suggest linear 
curves, even if there are no real power laws. 

It would thus be extremely useful if these empirical ob- 
servations could be backed by theoretical models where 
power laws can be either proven exactly or at least veri- 
fied beyond doubt by high statistics simulations. Unfor- 
tunately, equilibrium systems in general show anomalous 
scaling laws only at critical points which are codimcnsion 
one phenomena: One has to fine tune some parameter 
(e.g. temperature) to reach them, otherwise no power 
laws are obtained. Thus, they cannot be used to justify 
why such power laws should be seen in nature. 

A possible solution of this puzzle was indicated in |^|, |f| 
where it was suggested that many non- equilibrium sys- 
tems could be driven by their dynamics into a critical 
state. The main ingredients of this self organized criti- 
cality (SOC) is slow driving towards some instability and 
a mechanism to relax the tensions built up by the drive 
locally and partially. Since the tensions are not relaxed 
completely (in many models they are just redistributed), 
the state becomes marginally stable and apt to relaxation 
events on wider and wider length scales, which then lead 
to anomalous scaling laws. The paradigmatic model is 
the "sand pile" of j5j which does not describe real sand 
piles but which was proven exactly to show anomalous 
scaling laws at least for some of its observables || ^, || . 

Another model which was proposed to show SOC is the 
forest fire model introduced independently by Henley |J 
and by Drossel and Schwabl |l(|. In this lattice model 
with discrete time each site can be either occupied by a 
tree or empty. New trees are grown with small fixed rate 
p on empty sites, and with a rate / C p sites are hit 



by lightning strokes which then burn the entire cluster of 
trees connected to this site. Burning happens infinitely 
fast, so the only relevant parameter is the ratio 9 = p/f 
which also sets the scale for the average fire size (the 
number of trees burnt after one lightning) . Criticality is 
observed in the limit 9 — ► oo. 

The Drossel-Schwabl model (called DS model in the 
following) is different from other SOC models in two 
ways. 

• It involves not only the separation between two 
time scales (the slow build-up of stress and the fast 
relaxation), but involves three time scales: The fast 
burning of connected clusters of trees, the slow re- 
growth, and the even slower rate of lightnings. 

• The growth of trees does not lead to a state which 
is inherently unstable (as does the addition of sand 
grains to the top of a pile does), but only to a state 
susceptible to being burnt. Without lightning, the 
tree density in any patch of forest can go far beyond 
criticality. When the lightning strikes finally, the 
surviving trees have a density far above critical. 

Indeed, it was observed in |ll| that the stationary 
steady state of the model is not "critical" in most re- 
gions, in the sense that its tree density is not marginal 
for the spreading of fire. Rather, it is composed of large 
patches of roughly uniform density, most of which are ei- 
ther far below or far above the critical density for spread- 
ing. Nevertheless, power laws were observed for several 
observables, partly because these patches occur with all 
sizes, so that also the fires had a broad spectrum of sizes. 

While normal scaling with mean field exponents had 
been seen in [I0|, all subsequent simulations j^, 
|l3| , [l4| showed clear signs of anomalous scaling: 



11 



1- 



• The fire size distribution scaled, for small s (s is 
the number of trees burnt in one fire) and in the 



limit 9 — > oo, as P(s) ~ s 



with 



2.15; 



• For finite 9, P(s) is essentially cut off at s n 
with A ~ 1.1; 
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• The average rms. radius R(9) = (i? 2 ) 1 / 2 of all fires 
scaled as 9" with v ps 0.58. Here, R 2 is the Eu- 
clidean distance of a burning tree from the site of 
lightning, and the average is taken over all trees in 
all fires at fixed 9; and 

• The rms. radius R(s, 9) of fires of fixed size s scaled 
as s 1 / where the fractal dimension D of fires is 
D re 1.96 H to 2.0 [0. 

Finally, the average density of trees was found to be p = 
0.408- const/ V&. 

There were however seen already at that time large 
corrections to this scaling law and deviations from con- 
ventional ansatzes. Thus, 

• The determination of r was subject to large sys- 
tematic uncertainties |QJ ; 

• The scaling R ~ s 1 / 13 was observed in [[tlj only for 
fires of typical size (s ~ 9). For large s, R was 
significantly larger; 

• For finite 9, the fire distribution P(s) did not follow 
the normal scaling ansatz P = s 1_T </'(s/s max ) | pT| ; 

• There are (at least) two differently divergent length 
scales HJ]: The correlation length evaluated from 
all pairs of sites scales differently with 9 than R(9) 

©; 

• Finite size behaviour is abnormal |lq ]. 

In two recent publications, these problems were taken 
up again. In |17j it was claimed that they are due to 
non- leading corrections to scaling. In pl| a more radical 
solution was proposed with two distinct classes of fires 
which both contribute to the scaling limit. In the lat- 
ter paper also a connection to ordinary percolation was 
proposed, and the conclusion was backed by a "coarse 
grained" model which supposedly could mimic the DS 
model at extremely large 9 and on extremely large lat- 
tices. 

It is the purpose of this paper to present very large 
simulations which show quite unambiguously that none 
of these describe really the true critical behaviour of the 
DS model. While our simulations agree perfectly with 
previous ones for the lattice sizes and 9 values used there, 
we shall see that the supposed scaling laws are just tran- 
sients. Even the present simulations do not reach the 
true asymptotic regime, but some suggestive features of 
the true asymptotics do emerge. 

We describe our simulations and their results in the 
next section. In the last section we draw our conclusions. 

II. THE SIMULATIONS 

Our simulations are straightforward, with a few sub- 
tleties. They follow Ref. [jll| in taking p — > by making 
no attempts to grow new trees while a cluster burns. As 



in pi] we made exactly 9 growth attempts on randomly 
chosen sites between two successive lightnings, instead of 
letting this number fluctuate. In the large L limit this 
should not make any difference. But in contrast to [Q| , 
where a depth first algorithm had been used, we used a 
breadth first algorithm to burn the cluster |19(| . 

In order to simulate very large lattices with minimal 
storage, we use multi-spin coding, i.e. we use only one bit 
to store the status of a site. In this way we could simulate 
lattices of size Lx L, L = 65536 on computers with 1 GB 
main memory. Notice that we do not need to store for 
every tree whether it is burning or not, since the burning 
trees are stored in a separate list which is updated at 
each time step. Boundary conditions were helical, i.e. 
sites are indexed by one scalar index i, with neighbours 
having indices i ± 1 and i ± L, and with i + L 2 = i. The 
largest previous simulations |L7| had used L — 19000. 

We were careful to discard sufficiently long transients 
(between 1.6 6 and 1.2 x 10 7 lightnings; this is up to one 
order of magnitude longer than those in |l4|, [l7]]) before 
taking statistics. This is needed since excessively large 
fires occur during these transients, and thus the tail of 
the distribution P(s) is heavily influenced by them. We 
believe that previous analyses were affected by this prob- 
lem which is easily overlooked since bulk properties (such 
as p or P(s) for typical s) show much faster convergence. 
The total number of fires in each run used for averag- 
ing was between 10 9 (for 9 < 250) and 9.3 x 10 6 (for 
9 = 256000; see Table 1). Previous authors Q went 
only to 9 = 32768 with 2.5 x 10 7 fires. Compared to 
that, our statistics is larger by roughly 1 order of magni- 
tude. All simulations were done on fast Alpha worksta- 
tions. The CPU times per run varied between 15h and 
4 weeks. As random number generator we used Ziff's 
four-tap generator with period 2 9689 — 1 p0| . 

We tested for finite size effects by making for the same 
9 several runs on lattices of different sizes. Previously, 
finite size effects had been studied systematically in [[l6| , 
but only for much smaller lattices (L < 2000; the au- 
thors of |l^] called their analysis a finite size analysis, but 
they actually made a conventional scaling analysis with- 
out checking the finite size behaviour). For 9 ~ 64000, 
e.g., we made runs with L = 2 14 ,2 15 , and 2 16 . We veri- 
fied that distributions like P(s), R(s,9), or P(t) (t is the 
burning time of a fire) were independent of L within the 
statistical accuracy, i.e. any systematic L-dependence 
was masked by statistical fluctuations. 

Systematic dependencies were seen only for averaged 
quantities like (s), p, or R(9). Indeed, p can be mea- 
sured either immediately before or immediately after a 
lightning. Since each lightning burns in average (s) trees 
(here (s) is averaged over all lightnings, whether they hit 
a tree or an empty site), we have ^before — Pafter = (s)/L 2 . 
We found that pbeforo decreased with L, while p a fter in- 
creased with it (see Table 1). More precisely, we have 

Paftcr(i)«pM-0.6( S )/L 2 . (1) 

Notice that (s) is given for L — > oo exactly by (s) = 
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FIG. 1: (a) Log-log plot of P(s) versus s, for fixed 6 ranging 
from 125 to 256000. Here and in Figs. 2 and 5 we used loga- 
rithmic binning with 2 per cent bin size, to suppress excessive 
statistical fluctuations. 

(b) Enlarged part of the same data, multiplied with s 1 ' 19 . 



(1 - p)9 |(J. For finite L, one has to replace just p by 
(Pbefore + Patter )/2, to leading order in l/L. This was 
verified in the simulations, but it tests just stationarity 
and the absence of gross mistakes. 

The rms. radius of fires depended on L and 9 in a more 
complicated way, see Table 1. The data were less clear 
in that case, but they could be fitted by 



R(0) L=oo - R(6) L m 0.004 9 3/2 /L. 



(2) 



Results for the fire size distributions P(s) are shown 
in Fig.l. We see the approximate power decay s 1 ~ T for 
s < s max with s m ax roughly proportional to 9. But we 
see also the strong deviations from this power law first 
observed in jy]. Due to these deviations, P(s) decreases 
with 9 in the scaling region (making the effective value of 
r increase with 9), but has a growing bump at s rs s max - 
In Q it was conjectured that the asymptotic value of t, 
estimated from the scaling region 1 <C s <C 9 in the limit 
9 — > oo, is r = 2.45. This was based on heuristics and 
on simulations for small 9 on very small lattices (L 
1 .')()(): . In order to test it, we plotted P(s) in Fig. lb after 




FIG. 2: (a) Log-log plot of Pmt(s) versus s, for fixed 6 ranging 
from 125 to 256000. 

(b) Enlarged part of the same data, multiplied with s 019 . 



multiplying it with a suitable power of s. It is seen that 
s P(s) becomes flat in a wider and wider region of 
s as 9 increases. For very small fires (s < 300) P(s) 
decreases faster than s -119 , but there is no indication 
that the slope increases with 9 for 9 > 8000. Based on 
this evidence we would thus conclude that t = 2.19±0.01, 
ruling thereby out the value of - provided, of course, 
that there is no change of behaviour as 9 becomes even 
larger. 

In |]llf it was conjectured that the bump near s ~ 
Smax is due to the cut-off. If the integrated distribution 
-Pint(s) = Xi S '=s were just a power multiplied by 

a sharp cut-off, its derivative P(s) would have a bump 
where the cut-off sets in. This bump would consist of 
those events which would have been in the tail which is 
cut off. If this were right and Pi n t(s) would indeed show 
normal scaling, we should expect the height of the bump 
in Fig. lb to be independent of 9. But this is obviously 
not the case. Instead, it increases with 9, suggesting a 
different scaling law s 1_r for the envelope of the curves 
in Fig.la, with r' = 2.111 ± 0.006. 

A similar conclusion is reached by looking directly at 
fmt(s). Indeed, log- log plots of Pint (a) versus s are much 
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FIG. 3: Log-log plot of effective r values, obtained by drawing 
lines tangent to both maxima in Fig. 2b. Notice that this is 
only denned when P m t(s) is not convex, i.e. only for 9 > 1000. 
Error bars are roughly of the size of the symbols. 
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FIG. 4: Log-linear plot of 6> ' 582 R(0) versus 6, where the 
power of 6 was chosen such that the data fall roughly on a 
horizontal line. Without the factor 6l 582 i the data would be 
hardly distinguishable from a straight line on a log-log plot. 



more straight in the scaling region (Fig. 2a). But again 
a blow-up after multiplication by a suitable power of s 
shows that this is misleading (Fig. 2b). Even Pi n t(s) is 
not convex for 8 > 1000, and it develops an increasingly 
sharp shoulder near s — s max as 8 — > oo. We can try to 
obtain an alternative estimate of r by fitting straight lines 
such that they touch both maxima in Fig. 2b. Results 
of this are shown in Fig. 3. In contrast to the previous 
estimate of t they would indicate that the effective r 
decreases with 9 and is clearly less than 2.19 for 9 — > oo 
(also the value 2.159 of |L4j seems hardly compatible with 
the extrapolation of Fig. 3 to 8 — > oo). But a convergence 
to 2.11 seems quite possible, as we would expect if the 
bumps in fig. lb have a width independent of 8 (which is 
not excluded by the data). 

This obviously means that these scaling violations are 
not corrections which disappear in the limit 8 — ► oo, as 
was claimed in |Tfl ]. It rather indicates that the conven- 
tional scaling picture, and in particular the usual ansatz 

P(s) = fl 1 - T 0(«/w) J (3) 

is basically wrong, as was already conjectured in O, p"8| . 

Previous analyses indicated that s max actually in- 
creased faster than 9, roughly as 

Smax ~ 9 X (4) 

with A « 1.08. We verified this qualitatively, but verified 
also the finding of jnj that s max and thus also A are not 
well defined since the sharpness of the cut-off increases 
with 8. This is already obvious from Fig. 2b, but it per- 
sists for larger values of s not shown in this plot. It can 
also be seen by making copies of Fig. 2a on transparencies 
and overlaying them. 

Better scaling than for P(s) was seen in @ for R(9). 
Our present data are fully compatible with those of 
PI |l3| , |l4| , |l5[ ], but involve much higher statistics and 
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FIG. 5: Log-linear plot of R(s,0)/i/s versus s, for selected 
values of 9. 



cover a wider range of 8. Thus it might not be too sur- 
prising that we now do not see perfect scaling any more. 
But the observed deviations from a pure power law (see 
Fig. 4) are much larger than expected from subasymptotic 
corrections They clearly show that the previously 

seen power law was spurious and does not describe the 
asymptotic behaviour. A power law fit through the last 
two points would give R{8) ~ ^°- 563 j but this is obviously 
not yet the asymptotic behaviour. Again, as for P(s) and 
-Pint(s), we cannot yet say what the correct asymptotic 
behaviour will be. In any case, the claim of [[l4| that 
there are two different diverging correlation lengths be- 
comes obsolete. 

An indication of the origin of all these puzzles comes 
from looking at R(s,9). This is defined analogously to 
R(9), except that the averaging is now done over all fires 
of fixed size s. We expect R(s, 8) ~ s 1 ^ if fires are 
fractal with dimension D. Previous analyses had given 
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FIG. 6: Plot of v e fi versus 1/8, where i/ e g is defined by 
R(s,9) ~ s" ett at the inflection points in Fig. 5. 

D ps 2 (or slightly less). Thus we do not plot R(s,9) 
versus s, but R{s,9)/ s 1 / 2 . From Fig. 5 we see that there 
are 3 distinct regions which by and large coincide with 
the three regions (left bump, flat middle part, and right 
peak) in Fig. lb. 

• Region I: These are very small clusters, of size 
s < 300. For them R(s, 9) increases faster than 
y/s, which might suggest D < 2. But we rather 
interpret this as a finite-cluster artifact. Indeed, 
even compact clusters (with D — 2) with a 
crumpled boundary will show an effective fractal 
dimension < 2. These small clusters arise from 
lightnings which hit either very small regions with 
supercritical tree density which burn off com- 
pletely, or regions of subcritical density on which 
the fires form subcritical percolation clusters. 
Both mechanisms give compact clusters with fuzzy 
boundary. In any case, region I becomes less and 
less important as 9 increases. 

• Region II: This is roughly equal to the scaling 
region in fig. lb. Here R(s) is nearly proportional 
to y/s, i.e. D is very close to 2 - but not quite. 
Also, we see a clear decrease of the minimal slope 
in Fig. 5 with 9, but it seems not to be sufficient 
to give D = 2 in the limit 9 — ► oo (see Fig. 6). 
According to Fig. 6 the critical exponent v defined 
by R(s, 9) ~ s u in region II seems to converge to 
0.505 for 9 — > oo, corresponding to D = 1.98. This 
is half way between the best previous estimates 
D = 1.96 jl3| and D = 2.0 fn}. It indicates that 
clusters in Region II are fractal, but more compact 
than critical percolation clusters (D = 1.89). 

• Region III: In the region of very large fires, cor- 
responding to the increasing bumps in Fig. lb, the 
apparent fractal dimension decreases again. Unfor- 



FIG. 7: Plot of p versus l/y/d. These values were extrapolated 
to L — *• oo using Eq.(l), but for the largest L the extrapolation 
shifted them by less than the error bars (which are indicated 
inside the squares). The dashed line is a fit of the small 8 
values (8 < 8000) to an ansatz p = a — b/V0. 



tunately, due to the strong curvatures of the curves 
in Fig. 5 in this region, we cannot quote any definite 
dimension value. But by plotting R{s,9)/s v with 
suitable values of v on log-log scales, we see that 
(i) the maximal values of v increase slowly with 
9 (except for the very largest 9 where statistics is 
poor); and (ii) for our largest 9 and s, we have 
v Ki 0.64. Thus the largest fires are definitely more 
fractal than critical percolation clusters! 

In terms of the scenario with roughly homogeneous 
patches with constant tree density |ll], [l3| fires in 
region II correspond to single patches of typical size 
which are hit by lightning just when their density has 
reached about the critical percolation threshold. If light- 
ning would always strike exactly at criticality, and if all 
trees in the patch would have burnt during the last fire 
so that all trees now are placed randomly, this would give 
D = 1.89. Fires in region II are more compact mainly 
because it will take some time until a lightning strikes by 
chance the region, after it has reached the critical density. 

The larger a patch is, the bigger will be its perime- 
ter and thus also the chance that a fire "spills over" and 
burns also a neighbouring patch, and from this also a 
next patch, etc. Since this will happen only along parts 
of the perimeter, the resulting fires will be rather fuzzy, 
with effective dimension < 2. We propose that such fires 
dominate in region III. Although they are rather few in 
number, they are very important since they burn large 
parts of the entire lattice, and they lead to rearrange- 
ments of the global pattern of patches. Notice that fires 
in region III burn only small parts of the entire region 
they cover, leaving behind more unburnt trees than fires 
of type II. Since they dominate more and more as 9 in- 
creases, one might suspect that this leads to an increase 
of p with 9. 
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FIG. 8: Rescaled average fire sizes for fixed local density. 



FIG. 9: Distribution of local densities. 



Indeed, a slight increase of p with 6 had been seen in all 
previous simulations. The best previous estimates were 
p = 0.4075 — const/y/9 JllJ (unfortunately, the constant 
multiplying # -1 / 2 was estimated wrongly in JTl| ] due to a 
simple mistake) and p = (0.4084 ± 0.0005) - const/6 0A7 
fhi] , [l7| . The results of our present simulations are shown 
in Fig. 7. For 9 < 10 4 we see a perfect agreement with pre- 
vious results, but for larger 9 there is dramatic disagree- 
ment: Our measured values are higher than predicted by 
extrapolation from small 9, by up to 100 standard devia- 
tions. It is not clear why this was missed in |l7[| . But we 
might mention that no data for 9 > 10000 are shown in 
Fig. 1 of jl7j , although the authors claimed to have made 
high statistics measurements up to 9 — 32768. 

When plotted against log 9, our values of p follow 
roughly a straight line for 9 > 4000. We should of course 
not take this as the asymptotic behaviour, since p can 
never increase beyond p c = 0.5927. . . which is the crit- 
ical value for site percolation in 2 dimensions. But we 
can use it to obtain a very crude estimate of the order of 
magnitude when p w p c should be reached. It is 9 « 10 40 . 

Although any prediction based on such a large extrap- 
olation should be taken with great care, we conjecture 
that indeed p converges to p c for 9 — > 00. The main rea- 
son is that we see no other plausible scenario compatible 
with our present numerics. It is not clear whether fires 
in this limit correspond to critical percolation clusters, 
since we must expect that weak correlations in the tree 
densities survive in this limit, sufficiently so to spoil any 
agreement with uncorrelated percolation on large scales. 

A last hint in favour of our claim that large fires asymp- 
totically are dominantly associated to regions of critical 
tree densities (i.e., p ~ p c ), even though they might not 
form critical percolation clusters, is obtained by study- 
ing the mean sizes of fires which started in regions of 
given tree density. If a lightning hits a tree at site i, we 
define the local tree density p\ oc at this site as the num- 
ber of other trees in the surrounding square of size 9x9 
divided by 80. This size is of course rather arbitrary, 
but we obtained qualitatively similar results for squares 



of sizes 7x7 and 5x5. In Fig. 8 we plot the average 
fire size s(pi oc ), divided by the overall mean (s), against 
p\ oc . For small 9 we see a monotonic increase which is 
easily understood: Since even the largest patches of uni- 
form density are not much larger than the square, large 
fires can only result from regions with high local den- 
sity. This is no longer true for large 9. There, patches 
with very high density are probably small, otherwise they 
would already have been burnt down earlier. At values 
of 9 reached in this work the largest s(p\ oc ) are still for 
large /?i oc , but a pronounced peak at p p c develops 
where s(p\ oc ) has a local maximum. Since also the num- 
ber of sites with p\ oc w p c is much larger than those with 
Pioc S> p c (see next paragraph), we see that it is fires 
in regions with critical density which play an increasing 
important role as 9 — > 00. In any case, the very strong 
dependence of s(pi oc ) on 9 shows that we are still far from 
the asymptotic region where we expect this dependence 
to disappear. 

Finally, we show in Fig. 9 the distribution of local den- 
sities itself. It depends rather weakly on 9 (less than 
s(pioc), at least), but the precise way it does depend on 
9 is rather surprising and not yet fully understood. First 
of all, it develops an increasingly sharp maximum which 
slowly shifts from p\ ac « 0.4 to p\ oc « 0.6 as 9 increases. 
This is to be expected after our previous observations. 
What is unexpected and hard to explain is a shoulder at 
Pioc ~ 0.8 which develops for the largest 9 values. It is 
not very large but statistically highly significant (it was 
seen in all runs with 9 > 32000). Presumably related to 
it is a shoulder at small p\ oc (w 0.1) which also increases 
with 9: If patches with density > 0.7 burn, they leave be- 
hind extremely strongly depleted patches. One possible 
reason why such high density patches can survive at all is 
that many large fires are fractal and leave behind small 
disconnected regions of fairly high but subcritical den- 
sity. These regions then are too small to have a chance 
to be hit by lightning until their density has grown far 
beyond p c . 
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III. DISCUSSION 

The simulations reported in this paper leave little 
doubt that all scenarios proposed so far for the Drossel- 
Schwabl forest fire model do not describe the asymptotic 
behaviour for 6 — > oo, where the model should show SOC 
according to the standard folklore. Indeed, we do not see 
much indications for any power laws in this model, as all 
proposed scaling laws seem to be just transient. There 
are a number of observables which do show straight lines 
on log- log plots (such as Fig. la in the central region of s 
or the envelope to Fig. la), but it seems more likely that 
also these are spurious. 

This situation is of course not altogether new. There 
are a number of models which were supposed to show 
anomalous scaling, but closer scrutiny proved otherwise. 
A good example is the Bak-Chcn-Tang forest fire model 
p2| which at least in d = 2 is not critical |23| . Other ex- 
amples include the Newman-Sneppen model pi| where 
one can prove semi-numerically that the observed power 
laws are transient |^7j and maybe even the "classical" 
abelian sandpile model in d = 2. While power laws for 
waves were proven rigorously in that model, it might well 
be that the observed deviations from finite size scaling 
p5| , |26| do not herald multifractality but just simply no 
scaling at all. One indication for the latter is the fact that 



some scaling laws show violations which do not seem do 
vanish for increasing system sizes Also, some other 
quantities in the sandpile model which involve superposi- 
tions of many waves depend qualitatively on the geometry 
of the lattice (square vs. strip) ^7|. For a system with 
true scaling one would not expect this. 

The situation becomes even worse when going to real 
life phenomena. It does not seem unlikely that many of 
the observed scaling laws are just artifacts or transients. 
Problems of spurious scaling in models which can be sim- 
ulated with very high precision such as the present one 
should be warnings that not every power law supposedly 
seen in nature is a real power law. 
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critically reading the manuscript. 

Note added: After finishing this paper, I was informed 
that similar conclusions had been reached in a recent 
paper p8| , where the authors studied P(s) with compa- 
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paper. Unfortunately, neither the spatial sizes of fires 
nor precise values of p(9) were measured in (2§j. I am 
indebted to Gunnar Pruessner for sending me this pa- 
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9 


log 2 L 


N 


N trans 


density 


R(9) 


256000 


16 


9.3 x 10 6 


2.0 x 10 6 


0.410667 ± 0.000036 


2635.138 ± 2.59 


128000 


16 


7.5 x 10 6 


3.0 x 10 6 


0.409321 ± 0.000042 


1789.675 ± 1.96 




15 


15.2 x 10 6 


1.6 x 10 6 


0.409231 ± 0.000043 


1787.024 ± 1.37 


64000 


16 


46.4 x 10 6 


6.0 x 10 6 


0.407908 ± 0.000013 


1205.893 ± 0.55 




15 


42.5 x 10 6 


3.0 x 10 6 


0.407904 ± 0.000037 


1203.097 ± 0.56 




14 


37.7 x 10 6 


1.0 x 10 6 


0.407832 ± 0.000038 


1202.836 ± 0.59 


32000 


15 


42.3 x 10 6 


6.0 x 10 6 


0.406623 ± 0.000038 


808.948 ± 0.37 




14 


78.3 x 10 fi 


3.0 x 10 6 


0.406633 ± 0.000037 


807.253 ± 0.27 


16000 


15 


106.8 x 10 6 


8.0 x 10 6 


0.405428 ± 0.000009 


539.719 ± 0.16 




14 


103.3 x 10 6 


5.0 x 10 6 


0.405401 ± 0.000017 


539.063 ± 0.16 




13 


73.2 x 10 6 


2.0 x 10 6 


0.405360 ± 0.000029 


538.464 ± 0.20 


8000 


14 


166.3 x 10 6 


6.0 x 10 6 


0.404188 ± 0.000008 


359.049 ± 0.084 




13 


110.9 x 10 6 


3.0 x 10 6 


0.404135 ± 0.000028 


358.703 ± 0.102 


4000 


14 


207.8 x 10 6 


8.0 x 10 6 


0.402690 ± 0.000007 


239.213 ± 0.050 




13 


394.2 x 10 6 


4.0 x 10 6 


0.402691 ± 0.000010 


239.207 ± 0.041 




12 


332.7 x 10 6 


2.0 x 10 6 


0.402595 ± 0.000011 


238.992 ± 0.044 


2000 


13 


366.4 x 10 6 


6.0 x 10 6 


0.400614 ± 0.000005 


159.511 ± 0.026 




12 


525.7 x 10 6 


2.0 x 10 6 


0.400577 ± 0.000012 


159.441 ± 0.021 


1000 


13 


319.6 x 10 6 


9.0 x 10 6 


0.397636 ± 0.000005 


106.647 ± 0.018 




12 


458.4 x 10 6 


6.0 x 10 6 


0.397624 ± 0.000006 


106.681 ± 0.015 




11 


655.8 x 10 6 


4.0 x 10 6 


0.397570 ± 0.000010 


106.603 ± 0.013 


500 


12 


490.8 x 10 6 


6.0 x 10 6 


0.393419 ± 0.000006 


71.414 ± 0.0097 




11 


664.9 x 10 6 


3.0 x 10 6 


0.393399 ± 0.000010 


71.420 ± 0.0083 


350 


12 


663.4 x 10 6 


8.0 x 10 6 


0.390667 ± 0.000004 


58.093 ± 0.0068 


250 


12 


1000.0 x 10 6 


9.0 x 10 6 


0.387643 ± 0.000003 


47.757 ± 0.0045 




11 


664.4 x 10 6 


6.0 x 10 6 


0.387626 ± 0.000007 


47.768 ± 0.0056 




10 


1000.0 x 10 6 


2.0 x 10 6 


0.387574 ± 0.000008 


47.732 ± 0.0045 


175 


11 


733.6 x 10 6 


6.0 x 10 6 


0.383909 ± 0.000006 


38.788 ± 0.0043 


125 


11 


922.7 x 10 6 


8.0 x 10 6 


0.379837 ± 0.000006 


31.847 ± 0.0032 




10 


1000.0 x 10 6 


3.0 x 10 6 


0.379804 ± 0.000008 


31.850 ± 0.0030 


88 


11 


1000.0 x 10 6 


12.0 x 10 6 


0.374935 ± 0.000004 


25.903 ± 0.0025 




10 


1000.0 x 10 6 


8.0 x 10 6 


0.374910 ± 0.000005 


25.906 ± 0.0025 



TABLE I: Statistics and main results: N is the number of lightnings used for averaging, N tra ns that of lightnings discarded 
during the transients. The density is measured after fires got extinct and before new trees are grown. Notice that the density 
has much larger errors in some runs than in others of comparable statistics. This results from the fact that there are important 
long-ranged negative autocorrelations in the density time series, and I had not written out the information needed to take them 
into account in these runs. 



